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Study of gravitational-radiation induced merging rates of relativistic binary stars (double neutron 
stars; neutron star + black hole; double black holes) shows that the first-generation gravitational 
wave interferometers with an rms-sensitivity of 10~ 21 at frequency 100 Hz can detect 10-700 black 
hole and only ~ 1 neutron star coalescences in a 1-year integration time in a wide range of stellar 
evolution parameters. It is notable that modern concepts of stellar evolution predict that the first 
detection of gravitational wave will independently discover black holes. 

PACS number(s): 04.30.Db; 04.80.Nn; 97.60.Jd; 97.60.Lf 

The final merging stage of binary relativistic star evo¬ 
lution containing two compact stars (neutron stars (NS) 
or black holes (BH)) that merge due to gravitational- 
radiation induced orbit decay on a time-scale shorter 
than the Hubble time (1.5 x 10 10 years) are among the 
primary targets for gravitational wave (GW) interfer¬ 
ometers currently under construction (LIGO, VIRGO, 

GEO-600) [|IJ. Reality of such events in the Universe is 
confirmed by binary pulsar observations ||. To date, 5 
binary NS are known in our Galaxy; 3 of them should 
coalesce within the Hubble time ||. Based on the prop¬ 
erties of these binary pulsars, one may estimate the total 
number of such systems and hence the merging rate of 
binary NS in the Galaxy, and then extrapolate this rate 
to the volume of the Universe where GW-detectors will 
be able to detect GW-signal at a given signal-to-noise 
level. Such estimations made in the last five years ||] 
yield a “realistic” galactic rate of double NS mergings of 
10” 5 — 1CV 6 per year. 

Theoretical estimate of double NS merging rate may be 
deduced from stellar evolution H and give ~ 10 4 NS- 
NS coalescences per year per Galaxy. It can be shown 
that these two estimates do not contradict to each other 
if one takes into account NS+NS binaries which contain 
no pulsar (ejecting NS) and thus are unobservable by 
traditional radioastronomical means |(J . 

The stellar evolution theory also predicts that neutron 
star - black hole (NS+BH) and black hole black hole 
(BH+BH) binaries should be present in the Galaxy. So 
far the indirect existence of BH in stellar binary sys¬ 
tems has been argumented by a high mass (> 3 Mq) of 
the unseen companion in 11 close X-ray binary systems 
(?]. No PSR+BH binary system has yet been discov¬ 
ered in spite of optimistic theoretical expectations || of 
~ 1 per 1000 single pulsars. Theoretical estimates of 
NS+BH/BH+BH merging rate in the Galaxy are 1-2 or¬ 
ders of magnitude lower than NS+NS merging rate. The 
difficulty in doing such estimates is that additional pa¬ 
rameters of BH formation appears: the critical mass M cr 
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beyond which the star collapses into BH and the mass of 
BH formed. 

From the point of view of GW detection, NS+BH 
and BH+BH binary systems are advantageous in having 
higher ’’chirp mass” (M = M([i/M) 3 / 5 , M = Mi + M 2 , 
/r = M 1 M 2 /M) that determines the amplitude of the di¬ 
mensionless metric strain produced by a binary system: 
hamp ~ Ad 5 / 3 / 2//3 /r, where / is the GW frequency and 
r is the distance to the binary. In fact, matched filtering 
technique of data analysis |J allows enhancement of the 
signal-to-noise ratio from coalescing binaries S/N oc yfn 
where n = f 2 /f is the number of cycles of the signal 
passing through the detector frequency band A/ ~ /, 
so one usually considers the characteristic strain ampli¬ 
tude h c = h amp yfn ~ A4 5 / 6 / _1 / 6 /r. This means that at 
a given S/N the detector will be sensitive to a distance 
r oc Ad 5 / 6 . For example, if we consider a NS+NS bi¬ 
nary consisting of two identical NS Mi = M 2 = 1.4M© 
(A 4 ns « 1.2M 0 ) and a BH+BH binary consisting of two 
identical BH Mi = M 2 = 10M 0 (A4 bb ~ 8.7M 0 ), the 
limiting distance r bb — 5.2 r ns , or the volume available 
to search for such BH+BH systems is 140 times as large 
as the volume in which NS+NS binaries can be detected. 
Assuming homogeneous matter distribution, the detec¬ 
tion rate in a 1-year integration time with a given signal- 
to-noise ratio is 


N bh /N ns = (R bh /R ns ) x ( M bh /M ns ) 15/6 

where R bb and R ns are galactic merging rates for BH and 
NS binaries, respectively. Clearly, to answer the question 
what type of binary systems will be more numerous dur¬ 
ing 1-year operation time of a GW-interferometer one 
should study in detail different types of binary mergings 
in the Galaxy. 

In this paper we present calculations of NS+NS, 
NS+BH and BH+BH galactic merging rates in a wide 
range of main parameters of stellar evolution. We ad¬ 
dress the question how many binary mergings can one 
expect to detect with the rms-sensitivity h rms ~ 10 -21 
at f c = 100 Hz corresponding to S/N = 1 for GEO- 
600 and S/N = 3 for LIGO/VIRGO-I 0. We find that 
the prospect to discover BH-events (binary BH merging 
or BH+NS merging) is at least 10 times better than for 
NS+NS mergings. Moreover, for currently popular high 
recoil velocities of young pulsars (~ 400 km/s; |^0|) only 
BH+BH binaries may be detected by the first-order laser 
interferometers. 

Qualitatively, a crude estimate of the ratio of BH- 
containing binary merging rate to NS binary merging 
rate may be done as follows. BH is thought to result 
from the core collapse of massive stars. Single mas¬ 
sive stars lose roughly half their mass through inten¬ 
sive stellar wind; when in a close binary system, the 
mass may be transferred onto the secondary compan¬ 
ion. An approximate relation between the initial mass 
of the star and its core is M core ~ O.IJW^ [JO]]. Let 
us assume the mass of BH progenitor just before the 
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collapse be 35Mg, which would correspond roughly to 
M ms ~ 60 Mq on the initial main sequence. On the 
other hand, any star with M ms > 10Mq evolves to form 
a NS. Using the Salpeter mass function for star forma¬ 
tion ( dN/dM fa 1 (M/M©) -235 star per year), we obtain 
that BH formation rate relates to NS formation rate as 
(60/10) -1 ' 35 « 0.09. Extrapolating this logic to binary 
BH/NS systems, we might expect Rbh/Rns ~ 1/10, to 
a half-order accuracy. Actually, the situation is compli¬ 
cated by several factors: an asymmetry of the supernova 
explosion which may act more efficiently in the case of 
NS formation; mass exchange between the components; 
distribution by mass ratio, etc. All these factors will be 
accounted for in our calculations. 

To perform evolutionary calculations, we employ 
Monte-Carlo method for binary stellar evolution studies 
developed by us over last ten years (the Scenario Machine 
code); we refer to |l^| for a detailed description of the 
method and evolutionary scenarios used. The basic idea 
is to calculate the evolution of a representative number 
of binary stars (typically 10 6 ) whose orbital and physical 
parameters are distributed corresponding to the observa¬ 
tional data when available (e.g. initial masses, semimajor 
axes, recoil velocities of newborn neutron stars), or are 
chosen according to some model laws (e.g. initial rota¬ 
tion periods and magnetic fields of neutron stars). The 
calculated number of mergers is then scaled to the total 
stellar mass of the Galaxy (taken 1O 11 M 0 ) and after di¬ 
viding by the age of the Galaxy (1.5 x 10 10 years) yields 
the galactic merging rate R. Note (see below) that when 
reducing to the total mass of the Galaxy we must take 
into account the fraction of stars entering binary systems. 
This fraction is at least 50% 

A BH is known to be fully described by three param¬ 
eters: its mass Mbh, angular momentum, and electric 
charge. For our purposes, however, only mass is impor¬ 
tant as it determines the orbital evolution when the BH 
resides in a close binary system. At present, stellar-mass 
BH are thought to result from the core collapse of high- 
mass stars with M > M * where M* is the pre-supernova 
mass (see Q for recent discussion). It seems reasonable 
to assume that the mass of BH is proportional to M*, i.e. 
Mbh = kbhM *, 0 < kbh < 1. This mass in turn is calcu¬ 
lated from modern stellar evolution theory. Analysis of 
observational data on BH-candidates in binary systems 
H shows that the most plausible BH-formation parame¬ 
ters are kbh = 0.25 — 0.5, M* = 25 — 50. Some restriction 
to these parameters are discussed below. 

When in binaries, another important parameter ap¬ 
pears: the additional velocity Wbh imparted to BH during 
the anisotropic collapse. Here several different possibil¬ 
ities are feasible: (1) Wbh is proportional to the mass 
ejected during the collapse M e j = M*( 1 — kbh)', (2) the 
momentum MbhWbh is proportional to the ejected mass 
M e j\ (3) MbhWbh is proportional to the ejected momen¬ 
tum. 

We shall assume a universal mechanism giving 
anisotropic velocity for both NS and BH. In the case (1) 
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we have: 


Wbh = w ns (l - kbh)(l ~ M ov /M*) 1 ( 1 ) 

where Mov = 2.5 M 0 is the Oppenheimer-Volkoff limit 
for NS mass. This law is chosen assuming boundary con¬ 
ditions Wbh = 0 at kbh — 1 (he. when the total mass 
of the collapsing star goes into a BH) and Wbh = w ns 
once Mbh = Mov • Clearly, the case (2) is similar to the 
case (1) but produces lower BH-velocities while the case 
(3), instead, assumes Wbh = w ns irrespective of BH mass 
and thus produces higher BH velocities. We found, how¬ 
ever, that all three cases yield BH+BH/BH+NS merging 
rates practically coinciding with each other at high kick 
velocities, so here we represent the results for the case 
(1) only. 

The three-dimensional NS kick velocity is assumed to 
be arbitrarily directed in space and to be distributed so 
as to reproduce the observed pulsars’ transverse (two- 
dimensional) velocities © (see © for more detail): 

fLL(x)dx oc £ 0 ' 19 (1 + x e ' 72 )~ 1 ^ 2 dx (2) 

where x = w/wo, u>o is the characteristic velocity; the 
observed Lyne-Lorimer 2D-distribution is obtained at 
w o = 400 krn/s. The mean 3-D kick velocity differ from 
wq by ~ 10%. Our analysis of modern evolutionary sce¬ 
narios |L2| shows that the most probable kick velocities 
he in the range 200-400 krn/s. 

In Fig. |l] we plot the relativistic compact binaries’ 
merging rates as a function of the mean kick velocity as¬ 
suming Lyne-Lorimer distribution (]2|). BH-formation pa¬ 
rameters were taken within the limits M* = 15 — 50Mq, 
kbh = 0.25. Different scenarios for stellar evolution were 
considered (see ?? for more detail about evolutionary 
scenarios). From Fig. p] we see that the theoretical ex¬ 
pectation for the NS+NS merging rate in a model spi¬ 
ral galaxy with typical mass of 0.5 x 10 11 Mq in bi¬ 
nary stars lie within the range from ~ 3 x 10~ 4 yr _1 
to ~ 3 x 10 -5 yr _1 , depending on the assumed mean 
kick velocity and the shape of its distribution. Note that 
for spherically symmetric collapse the calculated rates 
coincide well with those found by Tutukov & Yungelson 
H who used another method of calculations. For Lyne- 
Lorimer law with the mean value of 400 krn/s, we obtain 
Rns+ns ~ 5x 10 -5 yr _1 . BH+BH and BH+NS merging 
rates are typically 1-2 orders of magnitude smaller. 
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FIG. 1. The mer g in g rate of NS+NS, NS+BH, and 
BH+BH binaries for Lyne-Lorimer kick velocity distribu¬ 
tion as a function of wo assuming BH formation parameters 
M t = 15 — 50 M©, kbh = 0.25, for different scenarios of binary 
star evolution in a model spiral galaxy with the total stellar 
mass 10 11 M© and 50% baryon fraction in binary stars. 


Two details in Fig. |l] are worth noting: 1) the binding 
effect at small kick velocities and 2 ) the smaller effect of 
high kicks on the BH+BH/BH+NS rate. The first fact is 
qualitatively clear: a high kick leads to the system disrup¬ 
tion; however, if the system has survived the explosion, 
its orbit would have a periastron always smaller than in 
the case without kick. During the subsequent tidal cir¬ 
cularization a closer binary system will form which will 
spend less time prior to the merging. The binding effect 
at small recoil velocities is more pronounced in the case of 
binary BH. At higher kicks their merging rate decreases 
slower than that of binary NS due to higher masses of 
the components. 

Having found galactic merging rates R , we can 
calculate them within a volume accessible for GW- 
observations of each type of binaries (which from the 
point of view of GW-signal differ, to the first approxi¬ 
mation, in having different chirp masses). Here we may 
either use normalization to the IR-luminosity per cubic 
megaparsec (as in Phinney’s paper [Q), or scale the galac¬ 
tic rates per baryonic density in stars and then extrapo¬ 
late it to the volume desired. In the latter case we obtain 
the formula 

U = 0.0063 R(e d /0.5) (H b /0.0046) h 2 75 Mpc~ 3 

where fib is the baryon density (in units of critical den¬ 
sity to close the Universe), /175 = H 0 /75 km/s/Mpc 
is the present value of the Hubble constant, ed is the 
fraction of baryons contained in binary stars. Typi¬ 
cally 0.25 < ed < 0.75. This normalization coincides 
with Phinney’s one assuming edflb = 0.0046. According 
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to H], fib = 0.0015 in stars within galactic disks and 
fib = 0.003 in spherical bulges of spiral galaxies and in 
elliptical galaxies. Calculations of NS+NS merging rate 
in elliptical galaxies show that although a strong de¬ 
crease in binary merging rate occurs, it tends to a nearly 
constant value several times lower than the binary merg¬ 
ing rate in a spiral galaxy with constant star formation. 
So taking edflb = 0.5 x 0.0046 would not be too far off 
and we use it to evaluate binary merging rates in a given 
volume of the Universe. Note, however, that in view of 
unknown exact value of all absolute rates we obtain 
may be uncertain to within a factor of 2. 

Consider now whether is it possible to restrict BH for¬ 
mation parameters from the existing astronomical obser¬ 
vations. From known 11 BH candidates in binary systems 
the evolutionary status of Cyg X-l is mostly well under¬ 
stood (a massive X-ray binary containing a 10 M 0 BH 
and a 30 Mg OB-supergiant). The evolution of 6 X-ray 
Novae - BH candidates is much more controversial. Only 
one CYg X-l-likc source is observed in the Galaxy, so for 
the number of such system we can adopt the lower limit 
N {Cyg X-l) > 1. On the other hand, as the analysis of 
stellar evolution shows ||, 1 binary PSR with BH should 
exist per 1000 single radiopulsars in the Galaxy. So far no 
such systems have been discovered, i.e. there is an upper 
limit 1V(PSR+BH)/^V(BPSR) < 1/700. Binary systems 
like Cyg X-l and PSR+BH are evolutionary related ||] 
and hence put bounds on BH formation parameters as 
shown in Fig. ||: kbh > 0.5, M* > 18M 0 (or equiva¬ 
lently, M cr > 35M 0 ). 


Presupernova Mass, M 0 
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FIG. 2. The galactic number of Cyg X-l-like systems (up¬ 
per curves) and the fraction of binary PSR with BH among 
the total number of single PSR in the Galaxy (lower curves) 
as a function of BH formation parameters. The upper hori¬ 
zontal scale is in presupernova masses M*, the lower scale in 
the critical main sequence mass M cr (both in solar masses). 
The fraction of presupernova mass to collapse into BH kbh is 
along the vertical scale. The shaded region shows the most 
plasuible BH formation parameters. 

Fig. | shows the expected total (NS+NS, NS+BH, 
BH+BH) and NS+NS detection rates at S/N = 1 level 
in 1-year integration on a GW-detector with the initial 
laser interferometers rms-sensitivity h rms = 10 -21 . All 
possible combinations of BH formation parameters M„, 
kbh with Lyne-Lorimer kick velocity law with wq = 400 
km/s were calculated for different scenarios of stellar evo¬ 
lution. It is seen from Fig. ||that NS+NS detection rate 
is ~ 0.3 — 0.7 per year whereas BH+BH/BH+NS co¬ 
alescences are parameters-dependent and can be much 
more numerous. The filled ”Loch-Ness monster”-headed 
region corresponds to the “most realistic” BH-formation 
parameters (M* > 18, kbh > 0.5, wq = 200 — 400 km/s 
and the “low mass-loss” scenario for evolution of single 
stars (see [?] for more detail). Within this region we may 
expect from ~ 10 to ~ 700 events per year, mostly (more 
than 80%) BH+BH coalescences. 



0.0 0.2 0.4 0.6 0.8 1.0 


kbh 

FIG. 3. The total merging rate of NS+NS, NS+BH, and 
BH+BH binaries as would be detected by a laser interfer¬ 
ometer with h rrrls = 10~ 21 for Lyne-Lorimer kick velocity 
distribution with w o = 200 — 400 km/s and BH progeni¬ 
tor’s masses M* = 15 — 50 M@, for different scenarios of 
binary star evolution as a function of kbh■ NS+NS merg¬ 
ings are shown separately. In all cases BH+BH mergings 
contribute more than 80% to the total rate. The filled 
“Loch-Ness-monster-head”-like region corrspond to BH for¬ 
mation parameters M* > 18 M@ and kbh = 0.5. 
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For a coalescing binary, the total mass M = M\ + M 2 
and the chirp mass A4 determine the final merging fre¬ 
quency (through the 3d Kepler’s law) and the GW- 
waveform amplitude, respectively. These mass distri¬ 
butions for binary BH mergings in the scenario with 
M* = 35 M 0 , kbh = 0.3 and wo = 400 km/s are shown 
in Fig. [|. The distributions are normalized so as to give 
J{dN/dM)dM = 1. 



FIG. 4. The total mass and chirp mass distributions for 
merging BH+BH binaries for M, = 35 Mg, kbh = 0.3 and 
Lyne-Lorimer kick with wo = 400 km/s. 

The results presented in this paper were obtained on 
the base of astronomically well-grounded binary stellar 
system evolutionary scenario. We then may ask the ques¬ 
tion: do there exist some parameters of the modern evo¬ 
lutionary scenario with which less than 1 coalescing rela¬ 
tivistic binary should be observed with our detector in 1- 
year integration time? The answer is “yes” only assuming 
kbh <0.1. Then low-mass BH are formed which are sub¬ 
jected to the disruptive action of both the high mass-loss 
during the collapse and the presumed collapse anisotropy. 
These small kbh seem to be unrealistic considering ob¬ 
served masses of BH-candidates ~ 10 Mq [|tJ (then the 
progenitor’s mass would be about 100 M@ which corre¬ 
sponds to implausibly high mass of 160 Mq on the main 
sequence). 

We conclude that for a wide range of BH-formation pa¬ 
rameters (M* ~ 18 80M 0 , kbh > 0.5) and a high birth 
velocity of newborn collapsed stars (300-400 km/s) we 
always expect ~ 10 — 700 binary merging events (mostly 
BH+BH pairs) to have the signal-to-noise ratio equal 
unity during 1-year integration by the first-order GW- 
detectors. Irrespective of the absolute number of merg¬ 
ing events (which are subjected to at least a factor of 2 
uncertainty due to unknown fraction of baryons in binary 
stars), the relative number of detected BH mergings at 
any detector should typically be ~ 10 times higher than 
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NS+NS coalescences. 


The authors acknowledge Prof. Leonid Grishchuk, Ed 
van den Heuvel and Bernard Schutz for useful discus¬ 
sions. KAP thanks the staff of the Department of Physics 
and Astronomy of University of Wales (Cardiff) for hos¬ 
pitality. The work was partially supported by the grant 
of Russian Fund for Basic Research No 95-02-06053a. 


[1] A. Abramovici et al, Science 256, 281 (1992); B.F. 
Schutz, in Les Houches Astrophysical School on Gravi¬ 
tational Waves, edited by J.-A. Mark and J.-P. Lasota 
(Cambridge Univ. Press, Cambridge, England, 1996, in 
press) (preprint AEI-003 February 1996). 

[2] J.H.Taylor, Phil. Trans. Roy. Soc. London A341, 117 
(1992). 

[3] D.J. Nice, R.W. Sayer, and J.H. Taylor, Astrophys. J. 
466, L87. 

[4] E.S. Phinney, Astrophys. J. 380, L17 (1991); S.J. Curran 
and D.R. Lorimer, Mon. Not. Roy. Ast. Soc. 276, 347 
(1995); E.P.J. van den Heuvel and D.R. Lorimer, preprint 
(1996). 

[5] J.P.A. Clark, E.P.J. van den Heuvel, and W. Sutantyo, 
Astron. Astrophys. 72, 120 (1979); V.M. Lipunov, K.A. 
Postnov, and M.E. Prokhorov, Astron. Astrophys. 176, 
LI (1987); R. Narayan, T. Piran, and A. Shemi, Astro¬ 
phys. J. 397, L17 (1991); A.V. Tutukov and L.R. Yun- 
gelson, Mon. Not. Roy. Ast. Soc. 260, 675 (1993); V.M. 
Lipunov, S.N. Nazin, I.E. Panchenko, K.A. Postnov, and 
M.E. Prokhorov, Astron. Astrophys. 298, 677 (1995); 
W.W. Dalton and C.L. Sarazin, Astrophys. J. 440, 280 
(1995); S.F. Portegies Zwart and H.N. Spreeuw, Astron. 
Astrophys. 312, 670 (1996). 

[6] . V.M. Lipunov, K.A. Postnov, and M.E. Prokhorov, sub¬ 
mitted to Mon. Not. Roy. Astron. Soc. (1996). 

[7] A.M. Cherepashchuk, Uspekhi Fiz. Nauk 166, 809 
(1996). 

[8] V.M. Lipunov, K.A. Postnov, M.E. Prokhorov, and E.Yu. 
Osminkin, Astrophys. J. 423, L121 (1994). 

[9] K.S. Thorne, in 300 Years of Gravitation, edited by S.W. 
Hawking and W. Israel (Cambridge University Press, 
Cambridge, England, 1987). 

[10] A.G. Lyne and D.R. Lorimer, Nature 369, 127 (1994). 

[11] I. Iben, Jr. and A.V. Tututkov, Astrophys. J. Suppl. 
54, 335; E.P.J. van den Heuvel, in Interacting Binaries, 
edited by S.N. Shore, M. Livio and E.P.J. van den Heuvel 
(Springer, Berlin, Germany, 1994). 

[12] V.M. Lipunov, K.A. Postnov, and M.E. Prokhorov, As¬ 
tron. Astrophys. 310, 489 (1996); Astrophysics and 
Space Physics Review, edited by R.A. Sunyaev, 9, part 
4, p.l (Harwood Academic Publishers, 1996) 

[13] P.Bender et al, in Particle and Nuclear Astrophysics and 
Cosmology in the Next Millenium, edited by E.W. Kolb 
and R.D. Peccei (World Scientific Publ., Singapore, 1995) 

[14] M. Fukugita, C.J. Hogan, and P.J.E. Peebles, Nature 


9 



381, 489 (1996). 

[15] V.M. Lipunov, K.A. Postnov, M.E. Prokhorov, I.E. 
Panchenko, and H. Jorgensen, Astrophys. J. 454, 593 
(1995). 


10 



